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ABSTRACT 



\A , Context. Cosmic rays (CRs) can be studied through the galaxy-wide gamma-ray emission that they generate when propagating in the interstellar 
» I ■ medium. The comparison of the diffuse signals from different systems may inform us about the key parameters in CR acceleration and transport. 
C3 ' Aims. We aim to determine and compare the properties of the cosmic-ray-induced gamma-ray emission of several Local Group galaxies. 

Methods. We use 2 years of nearly continuous sky-survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-ray 
Space Telescope to search for gamma-ray emission from M31 and M33. We compare the results with those for the Large Magellanic Cloud, the 
Small Magellanic Cloud, the Milky Way, and the starburst galaxies M82 and NGC253. 

Results. We detect a gamma-ray signal at 5<r significance in the energy range 200 MeV - 20 GeV that is consistent with originating from M31. 
The integral photon flux above 100 MeV amounts to (9.1 ± 1 .9 stat ± l.Ogys) X 10~ 9 ph cnT 2 s~'. We find no evidence for emission from M33 and 
derive an upper limit on the photon flux > 100 MeV of 5.1 X 10" 9 ph cm s~' (2cr). Comparing these results to the properties of other Local Group 
galaxies, we find indications of a correlation between star formation rate and gamma-ray luminosity that also holds for the starburst galaxies. 
Conclusions. The gamma-ray luminosity of M31 is about half that of the Milky Way, which implies that the ratio between the average CR densities 
in M31 and the Milky Way amounts to <f = 0.35 ± 0.25. The observed correlation between gamma-ray luminosity and star formation rate suggests 
that the flux of M33 is not far below the current upper limit from the LAT observations. 

Key words. Cosmic rays - Local Group - Galaxies: M31, M33, Milky Way, LMC, SMC, M82, NGC253 - Gamma rays: galaxies 



1. Introduction 

Cosmic rays (CRs) produce high-energy gamma rays through in- 
teractions with interstellar matter and radiation fields. The result- 
ing diffuse emissions directly probe CR spectra and intensities in 
galactic environments (e.g. Strong et al. 2007). The detection of 



the Small Magellanic Cloud (SMC; Abdo et al. 2010a) and de- 
tailed studies of the Large Magellanic Cloud (LMC; Abdo et al. 
2010b) and the Milky Way (MW; e.g. Abdo et al. 2009a) with the 
data collected by the Large Area Telescope (LAT) onboard the 
Fermi Gamma-ray Space Telescope enable comparative studies 
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of cosmic rays in environments that differ in star formation rate 
(SFR), gas content, radiation fields, size, and metallicities. 

Other galaxies in the Local Group that have been pre- 
dicted to be detectable high-energy gamma-ray emitters are M3 1 
(Andromeda) and M33 (Triangulum) due to their relatively high 
masses and proximity. So far, neither of these galaxies has been 
convincingly detected in high-energy gamma rays. M31 was 
observed by SAS-2 (Fichtel et al. 1975), COS-B (Pollock et 
al. 1981), and EGRET (Sreekumar et al. 1994), with the most 
stringent upper limit provided by EGRET being 4.9 x 10~ 8 
ph crrT 2 s _1 at a 95% confidence level (see Fig. 3 of Hartman 
et al. 1999). M33 has also been observed by COS-B (Pollock et 
al. 1981) and EGRET, providing an upper limit of 3.6 x 1(T 8 
ph crrT 2 s _1 (see Fig. 3 of Hartman et al. 1999). 

By comparing M31 properties to those of the MW, Ozel & 
Berkhuijsen (1987) estimated that the ratio £ of the CR density 
in M31 and in the MW is f ^ 1 and computed an expected 
>100 MeV flux fromM31 of 2.4 x l(T 8 £ph cirT 2 s -1 . Pavlidou 
& Fields (2001) made a comparable prediction of 1 x 10~ 8 
ph crrT 2 s _1 , based on the assumption that £ « 0.5, which they 
derived by comparing the estimated supernova rate in M3 1 and 
in the MW. Using the same approach, they also estimated the 
> 100 MeV flux of M3 3 to be 1.1 x 10~ 9 ph cirT 2 s _1 . 

If these estimates are correct, M3 1 should be detectable by 
the LAT after 2 years of sky survey observations, while M33 
still may fall below the current sensitivity limit. In this letter we 
report our searches for gamma-ray emission from M3 1 and M33 
with the LAT using almost 2 years of survey data. While we 
detect for the first time M31 just above the current sensitivity 
limit, we could only derive an upper limit for the flux from M33. 

2. Observations and analysis 

2.1. Data selection and analysis methods 

The data used in this work have been acquired by the LAT be- 
tween 8 August 2008 and 30 July 2010, a period of 721 days 
during which the LAT scanned the sky nearly continuously. 
Events satisfying the standard low-background event selection 
('Diffuse' events; Atwood et al. 2009) and coming from zenith 
angles < 105° (to greatly reduce the contribution by Earth albedo 
gamma rays) were used. Furthermore, we selected only events 
where the satellite rocking angle was less than 52°. We further 
restricted the analysis to photon energies above 200 MeV; below 
this energy, the effective area in the 'Diffuse class' is relatively 
small and strongly dependent on energy. All analysis was per- 
formed using the LAT Science Tools package, which is avail- 
able from the Fermi Science Support Center. Maximum likeli- 
hood analysis has been performed in binned mode using the tool 
gtlike. We used P6_V3 post-launch instrument response func- 
tions that take into account pile-up and accidental coincidence 
effects in the detector subsystems. 

2.2. M31 

For the analysis of M31 we selected all events within a rect- 
angular region-of-interest (ROI) of size 10° x 10° centred 
on (aj2000, feooo) = (00 h 42 m 44 s , +41°16'09") and aligned in 
Galactic coordinates. The gamma-ray background was modelled 
as a combination of diffuse model components and 4 significant 
point sources 1 that we found within the ROI. Galactic diffuse 

1 1FGL J0102.2+4223, 1FGL J0105.7+3930, 1FGL J0023.0+4453 
(Abdo et al. 2010c), and a hard source (T ~ 1.7) located at 
(a J2 ooo,<W) = (00 h 39 m 16\+43°27'07"). 



emission was modelled using an LAT collaboration internal up- 
date of the model gll_iem_v02 (e.g. Abdo et al. 2010c) refined 
by using an analysis of 21 months of LAT data and improved 
gas template maps with increased spatial resolution. Particular 
care was taken to remove any contribution from M31 and M33 
in the templates by excluding all gas with Vlsr < -50 km s _1 
within 2° x 3° wide boxes around (l,b) = (121°, -21.5°) and 
(l,b) = (133.5°, -31.5°) for M31 and M33, respectively 2 . In 
contrast to gll_iera_vQ2 3 , we did not include an E(B-V) tem- 
plate in the model because it includes some signal from these 
galaxies. We verified that the omission of the E(B-V) template 
did not affect the global fit quality over the ROI. The overall 
normalization of the Galactic diffuse emission has been left as 
a free parameter in the analysis. The extragalactic and residual 
instrumental backgrounds were combined into a single compo- 
nent assumed to have an isotropic distribution and a power-law 
spectrum with free normalization and free spectral index. The 
spectra of the 4 point sources were also modelled using power 
laws with free normalizations and free spectral indices. 

Figure 1 shows LAT counts maps for the energy range 200 
MeV - 20 GeV before (left panel) and after (right panel) sub- 
tracting the background model. For the purpose of highlighting 
emission features on the angular scale of M31, the counts maps 
were smoothed using a 2D Gaussian kernel of cr — 0.5°. In this 
representation, an elongated feature that roughly follows the out- 
line of M3 1 (as indicated by black contours) is already visible 
in the counts map before background subtraction (left panel). 
After this subtraction (right panel), this feature becomes the most 
prominent source of gamma-ray emission in the field. The re- 
mainder of the structure in the 'background subtracted' map is 
consistent with statistical fluctuations of the diffuse background 
emission, which illustrates that the signal from M3 1 is close to 
the actual detection sensitivity of the LAT. 

To test whether the emission feature is positionally consis- 
tent with M31, we performed maximum likelihood ratio tests 
for a grid of source positions centred on the galaxy. While the 
maximum likelihood ratio (or the maximum Test Statistic value 
TS; cf. Mattox et al. 2006) over the grid indicates the best- 
fitting source location, the decrease in TS from the maximum 
defines uncertainty contours that enclose the true source po- 
sition at a given confidence level. As usual, TS is defined as 
twice the difference between the log-likelihood of two alterna- 
tive models Xi and Xo, i.e. TS = 2(Xi - -£o)- Using a point 
source with a power-law spectrum, we obtain a best-fitting lo- 
cation of (aj2ooo.<*J20oo) = (00 h 42.4 m + 1.4 m ,+41°10' + 11') 
for the gamma-ray source, which encloses the centre of M31 
within the lcr confidence contour (quoted location uncertainties 
are at 95% confidence). Using instead of the point source an 
elliptically shaped uniform intensity region with a semi-major 
axis of 1.2°, a semi-minor axis of 0.3° and a position angle 
of 38° to approximate the extent and orientation of the galaxy 
on the sky 4 , we find a best-fitting location of (orj2ooo> feooo) = 
(00 h 43.9 m + 1.8 m , +41°23' + 22') that again encloses the centre 
of M31 within the lcr confidence contour. 

We determined the statistical significance of the detection, 
as well as its spectral parameters, by fitting a spatial template 
for M3 1 to the data on top of the gamma-ray background model 



2 For M31, the velocity cut left some residual in the Hi template 
owing to overlap in velocity with the MW along one side of M31. 

3 See the Galactic diffuse model description at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 

4 We estimated these parameters by adjusting an ellipse to the IRIS 
100 fim map of M31 (Miville-Deschenes & Lagache 2005). 
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Fig. 1. Gaussian kernel (cr = 0.5°) smoothed counts maps of the region of interest (ROI) in a true local projection before (left) and after subtraction 
of the background model (right) for the energy range 200 MeV - 20 GeV and for a pixel size of 0.05° x 0.05°. Overlaid are IRIS 100 jim contours 
of M3 1 convolved with the LAT point spread function to indicate the extent and shape of the galaxy. The boxes show the locations of the 4 point 
sources that have been included in the background model. 



that we introduced above. The M31 template was derived from 
the Improved Reprocessing of the IRAS Survey (IRIS) 100 /mi 
far infrared map (Miville-Deschenes & Lagache 2005). Far in- 
frared emission can be taken as a first-order approximation of 
the expected distribution of gamma-ray emission from a galaxy 
since it traces interstellar gas convolved with the recent mas- 
sive star formation activity. The spatial distributions of diffuse 
gamma-ray emission from our own Galaxy or the LMC are in- 
deed traced by far-infrared emission to the first order. From the 
IRIS 100 /mi map, we removed any pedestal emission, which 
we estimated from an annulus around M3 1 , and we clipped the 
image beyond a radius of 1.6°. 

Using this IRIS 100 //m spatial template for M31 and as- 
suming a power-law spectral shape led to a detection above the 
background at TS = 28.8, which corresponds to a detection sig- 
nificance of 5. Oct for 2 free parameters. We obtained a > 100 
MeV photon flux of (11.0 ± 4.7 stat ± 2.0 sys ) x 10~ 9 ph cirT 2 s" 1 
and a spectral index of F = 2 . 1 ± 0.2 stat + 0. 1 sys using this model. 
Systematic errors include uncertainties in our knowledge of the 
effective area of the LAT and uncertainties in the modelling of 
diffuse Galactic gamma-ray emission. As an alternative we fitted 
the data using the IRIS 60 /mi, IRIS 25 /mi, a template based on 
Ha emission (Finkbeiner 2003) or the geometrical ellipse shape 
we used earlier for source localization. All these templates pro- 
vide results that are close to (and consistent with) those obtained 
using the IRIS 100 /mi map. Fitting the data using a point source 
at the centre of M31 provided a slightly smaller TS (25.5) and a 
steeper spectral index (F = 2.5 + 0.2 stat + 0.1 sys ), which provides 
marginal evidence (at the 1 .8<x confidence level) of a spatial ex- 
tension of the source beyond the energy-dependent LAT point 
spread function. 

Using the gamma-ray luminosity spectrum determined 
from a GALPROP model of the MW that was scaled to the 



assumed distance of 780 kpc of M31 (Strong et al. 20 10) 5 
instead of a power law allows determination of the >100 MeV 
luminosity ratio r y between M31 and the MW. We obtain 
r y = 0.55 + 0.11 s tat ± 0.10 sys where we linearly added uncer- 
tainties in the assumed halo size of the model to the systematic 
errors in the measurement. The luminosity of M3 1 is thus about 
half that of the MW The model gives TS = 28.9, which is 
comparable to the value obtained using a power law, yet now 
with only one free parameter, the detection significance rises to 
5.3<t. According to this model, the >100 MeV photon flux of 
M31 is (9.1 ± 1.9 sta t + 1.0 sys ) x 10~ 9 ph cnr 2 s" 1 . 

We determined the spectrum of the gamma-ray emission 
from M3 1 independently of any assumption about the spectral 
shape by fitting the IRIS 100 /mi template in five logarithmi- 
cally spaced energy bins covering the energy range 200 MeV - 
50 GeV to the data. Figure 2 shows the resulting spectrum on 
which we superimposed the GALPROP model of the MW for 
r y - 0.55. Overall, the agreement between the observed spec- 
trum of M3 1 and the model is very satisfactory. The upturn in the 
spectrum at high energies, though not significant, could possibly 
be attributed to emission from the BL Lac object 1ES 0037+405, 
the only known blazar in the line of sight towards M3 1 . In a dedi- 
cated analysis above 5 GeV, we found a cluster of 6-7 counts that 
are positionally consistent with coming from that blazar. Adding 
1ES 0037+405 as a point source to our model and extending 
the energy range for the fit to 200 MeV - 300 GeV results in a 



5 We use throughout this work a representative model of the MW 
from Strong et al. (2010) with a halo size of 4 kpc and that assumes 
diffusive reacceleration. The model is based on cosmic-ray, Fermi-LKT 
and other data, and includes interstellar pion-decay, inverse Compton 
and Bremsstrahlung. Varying the halo size between 2 and 10 kpc affects 
the >100 MeV luminosity and photon flux by less than 10% and 3%, 
respectively. 
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Table 1. Properties and gamma-ray characteristics of Local Group and nearby starburst galaxies (see text). 



Galaxy 


d 


M m 




SFR 




L y 


q y 




kpc 


10 s M 


10 s M 


M Q vr" 1 


\Q- % ph cm" 2 s- 1 


10 41 phs" 1 


10~ 25 phs 'H-atom- 1 


MW 




35 ± 4< 7 > 


14 ± 2< 7 > 


1 - 3 <19) 




11.8±3.4< 28 > 


2.0 + 0.6 


M31 


780 ± 33 (1) 


73 ± 22 <8) 


3.6 ± 1.8 (14) 


0.35 - 1 (19 > 


0.9 + 0.2 


6.6+1.4 


0.7 + 0.3 


M33 


847 ± 60 <2) 


19 ± 8 (9) 


3.3 ± 0.4 <9) 


0.26 - 0.7< 20 > 


< 0.5 


< 5.0 


< 2.9 


LMC 


50 ± 2 {3) 


4.8 ± 0.2 (10) 


0.5±0.1 (15) 


0.20 - 0.25 <21) 


26.3 ± 2.0 (25) 


0.78 ± 0.08 


1.2 + 0.1 


SMC 


61 ± 3 (4) 


4.2 ± 0.4 (11) 


0.25±0.15 <16) 


0.04 - 0.08 (22) 


3.7 ± 0.7 (26) 


0.16 ±0.04 


0.31+0.07 


M82 


3630 ± 340 (5) 


8.8 ± 2.9 (12) 


5 ± 4 (17) 


13 - 33 (23) 


1.6 ± 0.5 <27) 


252 + 91 


158 + 75 


NGC253 


3940 ± 370 (6) 


64 ± 14< 13 ' 


40 ± 8 (18) 


3.5 - 10.4< 24) 


0.6 ± 0.4< 27 > 


112 + 78 


9 + 6 



References. (1) Stanek & Gamavich (1998); (2) Galleti et al. (2004); (3) Pietrzynski et al. (2009); (4) Hilditch et al. (2005); (5) Karachentsev et 
al. (2002); (6) Karachentsev et al. (2003); (7) Paladini et al. (2007); (8) Braun et al. (2009); (9) Grader et al. (2010); (10) Staveley-Smith et al. 
(2003); (11) Stanimirovic et al. (1999); (12) Chynoweth et al. (2008); (13) Combes et al. (1977); (14) Nieten et al. (2006); (15) Fukui et al. (2008); 
(16) Leroy et al. (2007); (17) Mao et al. (2000); (18) Houghton et al. (1997); (19) Yin et al. (2009); (20) Gardan et al. (2007); (21) Hughes et 
al. (2007); (22) Wilke et al. (2004); (23) Forster Schreiber et al. (2003); (24) Lenc & Tingay (2006); (25) Abdo et al. (2010b); (26) Abdo et al. 
(2010a); (27) Abdo et al. (2010d); (28) Strong et al. (2010): range based on GALPROP models with various halo sizes. 




1000 10000 
Energy (MeV) 



Fig. 2. Spectrum of the M31 emission obtained using the IRIS 100 jim 
spatial template. Red error bars are statistical, black error bars are sys- 
tematic uncertainties. The solid line shows an MW gamma-ray lumi- 
nosity model scaled to M3 1 and the dashed one a possible contribution 
of 1ES 0037+405 (see text). 



TS= 16-20 for the source, where the range reflects uncertainties 
in modelling the spectrum of the isotropic background compo- 
nent at energies >100 GeV. The fit suggests a hard power-law 
spectral index (F = 1.2 + 0.4), which explains why the source 
is only seen at high energies. Within 200 MeV - 20 GeV, how- 
ever, the source contributes only ~8 counts, a number that is tiny 
compared to the ~240 counts that are attributed to M3 1 . The im- 
pact of 1ES 0037+405 on the flux and gamma-ray luminosity 
estimates for M3 1 is thus negligible. 

We also repeated our analysis for a larger ROI of size 
20° x 20° in which we found 14 point sources in our LAT in- 
ternal source list. Searching for the faint signal from M31 in 
such a large ROI relies on the accurate modelling of the spatial 
distribution of the diffuse gamma-ray background over a large 
area, which is an important potential source of systematic uncer- 
tainties. Nevertheless, results obtained for this large ROI were 
consistent with those obtained for the 10° x 10° ROI. 

2.3. M33 

For the analysis of M33 we selected all events within a rect- 
angular ROI of size 10° x 10° centred on (aj2ooo, feooo) = 



(01 h 33 m 51 s ,+30°39'37") and aligned in Galactic coordinates. 
Within this field we detected 3 background point sources 6 that 
we included in the background model. The remainder of the 
analysis was similar to what was done for M3 1 . 

We did not detect any significant signal towards the direction 
of M33. Using a spatial template based on the IRIS 100 fim map 
of M33 and taking the GALPROP models of Strong et al. (2010) 
for the spectral shape, we derived an upper > 100 MeV flux limit 
of 5.1 x lO^phctrT 2 s- 1 (2<r). 

3. Discussion 

Based on the flux F y measured for M31 and the flux upper 
limit for M33, we computed the >100 MeV photon luminosi- 
ties L y = 4nd 2 F y and average emissivities q y = L y /Nu, which 
we compare to the values obtained for the MW, the LMC, and 
the SMC (see Table 1). Here, d is the distance of the galaxy and 
iV H = 1.19 x 10 57 (M H i + Mh 2 ) is the total number of hydrogen 
atoms in a galaxy, with Mm and Mh, in units of M . Quoted 
uncertainties in L y and q y include uncertainties in distance and 
hydrogen mass of the galaxies. The variations in L y and q y from 
one galaxy to another may inform us about how the CR popula- 
tion is affected by global galactic properties. From the q y values, 
we estimate the ratio £ of the average CR density in M3 1 and 
in the MW to £ = 0.35 + 0.25, consistent with the estimate of 
Pavlidou & Fields (2001). On the other hand, the flux upper limit 
for M33 allows for an average CR density in that galaxy that is 
above the MW value, hence up to a few times greater than the 
£ = 0.2 estimated by Pavlidou & Fields (2001). 

By comparing the L y of our sample of Local Group galax- 
ies to their total hydrogen masses and SFRs, we find a close 
correlation between L y and SFR and greater scatter between L y 
and gas mass (see Fig. 3). In the bottom panel of Fig. 3, the 
ranges of SFR values, which have been rescaled to the distances 
d adopted here, reflect uncertainties in the SFR estimates based 
on the various methods used to determine them (see Table 1). 
There is a clear trend toward increasing L y with increasing SFR, 
with L y = (7.4 + 1.6) x SFR L4±a3 when fitted by a power law, 
where L y and SFR are in units of 10 41 phs 1 and M yr _1 , re- 
spectively. We also added the luminosities derived by Abdo et al. 
(2010d) for M82 and NGC253 to this plot, illustrating that the 

6 1FGL J0134.4+2632, 1FGL J0144.6+2703, and 1FGL 
J01 12.9+3207 (Abdo et al. 2010c). 
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Fig. 3. Gamma-ray > 100 MeV luminosity versus total number of hy- 
drogen atoms (top panel) and star formation rate (bottom panel) for 
Local Group galaxies and the starbursts M82 and NGC253. In the bot- 
tom panel, the lines are power-law fits to the data for the MW, M31, the 
LMC, and the SMC, for which the slope was free (solid) or fixed to 1 
(dashed). 



relation obtained for Local Group galaxies also holds for nearby 
starburst galaxies. Assuming that it also holds for M33 allows 
estimation of the luminosity of Ly ~ (1 — 4) X 10 41 phs~' for 
this galaxy, corresponding to a > 100 MeV flux of (1 - 4) X 10~ 9 
ph ctrT 2 s _1 . M33 thus may be within reach of the LAT within 
the next few years. 

The L y -SFR plot does suggest a correlation in common for 
Local Group and starburst galaxies. Although it is premature 
to draw conclusions about any strong correlation over such a 
wide range of galaxy properties because of the small size of our 
sample, if such a correlation exists, it would be analogous to 
the well-known tight correlation between radio and far-infrared 
emission over a wide range of galaxy types (e.g. Murphy et 
al. 2006). The latter is linked to the relation between CRs and 
SFR, and although not yet fully understood, it is thought to re- 
sult to some extent from CR electron calorimetry. While pro- 
ton calorimetry clearly can be excluded as an explanation of 
the L r -SFR correlation because the intermediate-size galaxies 
of the Local Group are thought to be very inefficient at retaining 
CR protons, the dominant CR component (Strong et al. 2010), 
a correlation may relate to the contribution of CR leptons to 
the gamma-ray emission. Depending on the ISM and CR trans- 
port conditions, CR leptons may lose their energy predominantly 
through gamma-ray-emitting processes (like inverse-Compton 
or Bremsstrahlung, as opposed to ionization and synchrotron) 



and dominate the total gamma-ray luminosity 7 . This could drive 
the correlation between L y and SFR for galactic systems with 
high lepton calorimetric efficiency. Whatever the explanation for 
this global correlation, it is worthwhile noting that it holds de- 
spite the fact that conditions may vary considerably within a 
galaxy (e.g. the peculiar 30 Doradus region in the LMC, or the 
very active cores of starbursts). 

The L y vs SFR plane therefore seems to hold potential for 
defining constraints on CR production and transport processes. 
The inferred L y values are, however, not uniquely due to CR- 
ISM interactions but include a contribution of individual galactic 
sources such as pulsars and their nebulae. The relative contribu- 
tions of discrete sources and CR-ISM interactions to the total 
gamma-ray emission very likely vary with galaxy properties like 
SFR, which may complicate the interpretation of any L y trend in 
terms of CR large-scale population and transport. 

Also more exotic processes, such as annihilation or decay 
of WIMPs (weakly interacting massive particles), might con- 
tribute to the overall signal from M31. Several extensions of 
the Standard Model of particle physics naturally predict the ex- 
istence of WIMPs (e.g. supersymmetry, universal extra dimen- 
sions). Rather than focusing on a specific scenario, we estimate 
a conservative upper bound on this contribution in the case of 
a generic 100 GeV WIMP annihilating exclusively into bottom 
quarks, which is one of the leading tree level annihilation chan- 
nels of a WIMP predicted by supersymmetric theories. The nor- 
malization of the predicted spectrum is initially set to zero and is 
increased until it just meets, but does not exceed, the 95% confi- 
dence upper limit on the measured M3 1 spectrum at any energy. 
We find that when assuming an Einasto dark matter halo pro- 
file (Navarro et al. 2010) that matches the M31 kinematic data 
(Klypin et al. 2002), this contribution corresponds to a 95% con- 
fidence upper limit on the annihilation cross section of approxi- 
mately 5x10" 



1 25 cm 3 s 1 
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Appendix A: Gamma-ray spectrum of M31 

Table A. 1 provides the intensity values of the M3 1 gamma-ray 
spectrum that is shown in Fig. 2. Statistical errors are at the 
lcr confidence level, and the upper limit for the 16.6 - 50.0 
GeV energy bin at the 2<x confidence level. Systematic errors 
include uncertainties in our knowledge of the effective area of 
the LAT and uncertainties in the modelling of diffuse Galactic 
gamma-ray emission. The former were determined using mod- 
ifications of the instrument response functions that bracket the 
uncertainties in our knowledge of the LAT effective area. The lat- 
ter were determined by deriving spectra for variations of the dif- 
fuse Galactic models that make use of either an E(B-V) template 
or for which the gas templates have been replaced by the IRIS 
100 yum map, from which emission associated to M31 has been 
removed. Both types of systematic uncertainties were added lin- 
earily. 

The last column gives the number of counts attributed to 
M3 1 in each of the energy bins from the fit of a spatial model 
to the present data. 

Table A.l. Measured spectrum of M31 (see text). 



Energy Intensity Stat, error Sys. error Counts 
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0.15 


15.9 


5493 -16572 


1.27 


0.71 


0.26 


7.0 


16572 - 50000 


< 12.5 




3.6 


< 24.4 



